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Abstract 

We consider a cosmological scenario in which a scale-invariant spectrum of curvature per- 
turbations is generated by a rapidly-evolving equation of state on a slowly expanding back- 
ground. This scenario generalizes the "adiabatic ekpyrotic" mechanism proposed recently in 
arXiv:0910.2230. Whereas the original proposal assumed a slowly contracting background, 
the present work shows that the mechanism works equally well on an expanding background. 
This greatly expands the realm of broader cosmological scenarios in which this mechanism 
can be embedded. We present a phase space analysis and show that both the expanding and 
contracting versions of the scenario are dynamical attr actors, with the expanding branch 
having a broader basin of attraction. In both cases, a finite range of scale invariant modes 
can be generated within the regime of validity of perturbation theory. 



1 Introduction 



There is mounting observational evidence that the large scale structure originated from a nearly 
scale invariant and nearly Gaussian spectrum of primordial density perturbations. While these 
statistical properties are consistent with the simplest inflationary models [1], a critical question 
for early universe cosmology is whether inflation is unique in making these predictions. This has 
motivated the quest for alternative scenarios, from the pre-big bang scenario [2], to string gas 
cosmology j3] — [7j, to ekpyrotic theory [8j — [39J. (See [381 [39] f° r reviews of ekpyrotic cosmology.) 

A zeroth-order benchmark for a successful theory of the early universe is to explain the observed 
homogeneity, isotropy and flatness of our universe. There are only two known cosmological phases 
that make the universe increasingly homogeneous, isotropic and flat. The first is of course inflation, 
characterized by accelerated expansion after the big bang. This requires a scalar field (or fluid) 
with negative equation of state w < — 1/3. A second possibility is ekpyrosis, a phase of slow 
contraction before the big bang [181 Ell HO]- This corresponds to a scalar field with stiff equation 
of state w ^> 1. The smoothing power of ekpyrotic contraction was recently confirmed using 
numerical relativity simulations [UJ. 

Another benchmark is to generate a nearly scale invariant and Gaussian primordial spectrum. But 
even this is not sufficient. Inflation not only generates perturbations with the desired properties, 
but it does so within a cosmological background that is a dynamical attractor. Indeed, on super- 
horizon scales the curvature perturbation on uniform-density slices (12] — [44J, £ m 8a/ a, measures 
differences in the expansion history of distant Hubble patches [H]. Since £ — > const, at long 
wavelengths in (single-field) inflation, the background is an attractor |45j . 

Achieving both scale invariance and dynamical attraction in alternative scenarios has proven 
challenging. The mode function equation for ( in a contracting, matter-dominated universe takes 
an identical form as in inflation [45 1 147 1 li"8"] ; but ( grows outside the horizon in this case, indicating 
an instability. This is not surprising, since an anisotropic stress component, for instance, will 
blueshift faster than dust. Similarly, contracting mechanisms that rely on a time-dependent sound 
speed are inevitably unstable [49] . The contracting phase in the original ekpyrotic scenario [H], 
with V((j)) = —VQe~ c ^^ Mpi , is an attractor [TBI 124] . as mentioned above, and correspondingly the 
growing mode for ( is a constant in the long wavelength limit. However, the resulting spectrum 

1 



for ( is strongly blue [121 H31 HEJ EI], rather than scale invariant. A scale invariant spectrum 
can be obtained through entropy perturbations [151 EH 1251 [251 122], as in the New Ekpyrotic 
scenario [261 EH EH], but this requires two scalar fields. Even in this case, the entropy direction is 
generically tachyonically unstable [301 El] ■ 

The recently-proposed adiabatic ekpyrotic mechanism [351 EBl E3 [53] , based on a single scalar field 
with canonical kinetic term, offers a counterexample. The mechanism can be realized with fairly 
simple potentials, such as the "lifted exponential" 

V{<j>) = V Q (1 - e- c ^ M n , (1) 

with Vo > and c> 1. This potential takes the form of a plateau that has been lifted to positive 
energy at large positive <fi with a steep waterfall around = 0. See Fig. [TJ The key difference 
compared to earlier renditions of ekpyrotic cosmology is that the potential assumes positive values 
for part of the evolution. 

Scale invariant adiabatic perturbations are generated during a transient phase in which the equa- 
tion of state parameter, 

e = + (2) 

grows rapidly from e< 1, where the constant term dominates, to e ~ c 2 /2 3> 1, where the negative 
exponential term dominates. During this transition, the scale factor is nearly constant, while the 
equation of state parameter varies rapidly as e ~ 1/r 2 , where r is conformal time. This is referred 
to as the transition phase. The quantity z = a(r) -y/2e(r), which determines the evolution of £, 
therefore scales as z ~ (— t) 1 — exactly as in inflation, where e w const, and a(r) « 1/ (— r)! The 
two-point function is, therefore, identical to inflation. (Another counterexample proposed recently 
relies on a rapidly- varying, super luminal sound speed c s (r) [571 EH1 EH]. See [HH EH EH [62] for 
earlier related work.) 

More generally, it has recently been shown that in fact there are only three single-field cosmologies 
with unit sound speed capable of generating a scale invariant spectrum for ( on an attractor 
background [5U[55]: i) inflation, with a(r) ~ and e ~ constant; ii) the adiabatic ekpyrotic 
mechanism [351 EZ] mentioned above, with e ~ 1/r 2 on a slowly contracting background; and Hi) a 
novel version of the adiabatic ekpyrotic mechanism, in which the background first slowly expands, 
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Figure 1: Depiction of the "lifted exponential" potential, V((f>) = Vb(l — e c< ?V a/ pi). At large field 
values the potential is nearly constant, and there is a steep waterfall around = 0. 

then slowly contracts. The analysis has been generalized to the case of a time-dependent sound 
speed (S3 EH]- See [5TJ E21 E3] for related work. 

Although all three cosmologies yield identical two-point function, the degeneracy is broken by 
the three-point function. In contrast with the extremely Gaussian spectrum of the simplest infla- 
tionary models, the rapid variation of e in the adiabatic ekpyrotic mechanisms leads to strongly 
scale-dependent non-Gaussianities, which peak on small scales [37]. For the lifted exponential 
potential 0, this results in a breakdown of perturbation theory on small scales, both classically 
and quantum mechanically |37J. As shown in |37J, however, these pathologies, and related issues 
raised in [36J, can be avoided by considering more general forms of the potential, such as 

V (0) = V (1 - e -<^' Mw ) , (3) 

where the exponent decreases from c to a much smaller value b <C c after an acceptable range of 
scale invariant modes has been generated. The resulting spectrum is perturbative on all scales, 
but, because of the modified evolution, is now scale invariant over a finite range of modes, spanning 
a factor of 10 5 in k space, or a dozen e-folds. While limited, this range is sufficient to account for 
observations of the cosmic microwave background and the large scale structure. 



Earlier analyses of the adiabatic ekpyrotic mechanism have focused on the case where the universe 
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is contracting throughout. In this paper, we instead explore the phenomenology and implications 
of the expanding case, consisting of an "expanding transition phase", followed by a contracting 
ekpyrotic scaling phase. In doing so, we are motivated by the issue of embedding this mechanism 
in broader cosmological scenarios. Specifically, one attractive feature of the expanding branch is 
that we expect its basin of attraction to be much broader than in the original contracting version. 
Indeed, while the contracting transition phase in the original adiabatic ekpyrotic mechanism is a 
dynamical attractor, the evolution prior to the transition phase is generally not. At sufficiently 
early times, the field lies on a flat plateau of its potential, and, because the universe is contracting, 
any additional amount of kinetic energy is blueshifted and threatens to dominate the energy 
density. In the initially expanding version studied here, any additional kinetic energy will instead 
redshift away, and we therefore expect the solution to be stable for all times. 

After reviewing the contracting version of the adiabatic ekpyrotic mechanism in Sec. [2j we show 
in Sec. [3] that this mechanism works equally well on a slowly expanding background and gen- 
erates a scale invariant spectrum of curvature perturbations. Since scale invariance relies on a 
rapidly- varying equation of state parameter while the scale factor is nearly constant, the density 
perturbation spectrum is to a good approximation insensitive to whether the background is ex- 
panding or contracting during mode production. We check this explicitly in Sec. [4] by numerically 
integrating the mode function equation in both the expanding and contracting cases. The resulting 
spectra are indistinguishable (see Fig. [2] for a preview). Section [5] revisits the analytic calculation 
of the background evolution, allowing for more general initial conditions. We quantify the extent 
to which the scale invariant phase shortens as a function of initial scalar field kinetic energy. In 
Sec. [6j we study the issue of stability globally by performing a phase space analysis, including a 
wide range of initial conditions, both for the expanding and the contracting branch. This confirms 
that the transition phase is an attractor in both cases, as indicated at the perturbative level by £ 
having a constant growing mode as k — > 0. More broadly, this analysis also shows that the basin 
of attraction is broader in the expanding case, as anticipated in the previous paragraph. In Sec. [7] 
we study the connection to an inflationary precursor to the transition phase, which follows in the 
expanding case from trusting ([I]) at sufficiently large field values. We briefly summarize our main 
results and discuss prospects for future directions in Sec. [8j 

The expanding adiabatic ekpyrotic mechanism studied here suffers from the same strong coupling 
issues as the original contracting version. In particular, the calculation of the three-point amplitude 
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in [37] neglected the time-dependence of the scale factor, and hence should apply to the expanding 
case as well. (One technical difference in our case is that e momentarily blows up at the end of the 
transition phase, because H = at that time. As the numerical analysis of Sec. [4] clearly shows, 
however, the two-point function is insensitive to this momentary divergence. We expect that the 
same is true for the three-point function.) In particular, just like in the contracting branch, strong 
coupling can be avoided by considering more general forms of the potential, such as ([3]). We will 
not repeat the discussion of non-Gaussianities in this paper, and we refer the interested reader 
to [37] for more details. 



2 Review of the Adiabatic Ekpyrotic solution 

It is instructive to review the mechanism presented by [351 I3~7] in which scale-invariant curva- 
ture perturbations are generated by a fast-rolling scalar field during a slowly contracting phase. 
Consider a canonical scalar field coupled to gravity 

S = j d'x^g - \sTd^ - W)) , (4) 

where Mpj 2 = 8nG. For concreteness, we focus on the "lifted exponential" potential (fl|. This 
form of the potential is simplest for analytical calculations but is not meant to represent a realistic 
scenario. Indeed, to avoid strong coupling issues [36], realistic models must consider a more general 
potential of the form V ((f)) = Vq (l — e~ c ^^' Mpi ) , where c(<f>) decreases to a smaller value after 
a suitable range of scale invariant modes has been generated. See [37] for details. Since the 
production of these scale-invariant modes occurs as the field traverses a small range A(f> -C Mp\, 
there is a reasonable amount of freedom in specifying the global form of the potential. 

The evolution of the scalar field on a cosmological background is governed by the equation of 
motion 

+ 3H<f) + = . (5) 

During the phase of interest, the background evolves extremely slowly so we may ignore Hubble 
friction. As shown in [35], this approximation can be rigorously justified a posteriori by explicitly 
computing the first order correction and verifying that it is indeed small. Intuitively, the field 
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falls down a steep waterfall, hence the transition occurs within a Hubble time. This "fast roll" 
approximation leads to the simplified equation 



which is solved by 



where — oo < t < 0. In writing down this solution, we have assumed that the field initially has 
negligible kinetic energy at early times, hence the total energy is Vq. In Sec. [5j we will consider 
departures from this choice and the impact on the spectrum of perturbations. 

To solve for the evolution of the background, we substitute (|tJ) into H = —<fi 2 /2Mp l to obtain 
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H(t) = ^7 -H Q . (8) 



At sufficiently early times, the Hubble parameter is approximately constant, H(t) — Hq, where 
the integration constant Hq is fixed by the Friedmann equation, 

SM^H'-^ + V^^Vq. (9) 

The authors of [35J focused on the contracting branch, 



H ° = ^m- < 10 > 

Indeed, since Hq > and t is negative, H(t) is manifestly negative, corresponding to a contracting 
universe. (Note that the present sign convention differs from [35J, where H(t) = 2/c 2 t + Hq and 
Hq is a negative quantity.) 

The phase of interest is the regime where the Hq term dominates the Hubble parameter, H(t) ~ 
Hq, in which case e = —H/H 2 ~ 1/t 2 . Additionally, the assumption that the background remains 
nearly static implies a(t) ~ 1. The quantity z = a(r) -y/2e(r), which determines the evolution of 
C therefore satisfies 

z = aV2e ~ - , (11) 

T 

exactly as in inflation, where e ~ const, and a(r) ~ l/(— r). Moreover, as in inflation, the growing 
mode for ( goes to a constant, indicating that the background is a dynamical attractor [45J. The 
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two cosmologies yield identical power spectra, and therefore can be considered "dual" to one 
another at the level of the two-point function. 

The transition phase ends when H(t) ~ const, is no longer a good approximation. Subsequently, 
the Hubble parameter tends to an ekpyrotic scaling 

H ek (t) « A , (12) 

while the scale factor slowly decreases as 

0ek (t) ~ (-i) 2/c2 . (13) 

This ekpyrotic scaling phase is a necessary component of the story. Because the Hubble horizon is 
nearly static during the transition phase, the scale invariant modes created in this phase remain 
inside the Hubble radius. A subsequent phase is therefore necessary to push these modes outside 
the Hubble horizon. The scaling ekpyrotic phase fills this role — since the universe is slowly 
contracting, modes are gently pushed outside the Hubble horizon without disturbing their spec- 
trum. (Eventually, the ekpyrotic phase must itself terminate before the big crunch. We envision 
that it is followed by a bounce to an expanding, radiation-dominated phase. At the level of a 
four-dimensional effective theory, a stable non-singular bounce can be achieved either through a 
phase of ghost condensation [62] , as in the New Ekpyrotic scenario, or through a phase of galileon 
domination [64]. See [65] [66] for recent supersymmetric extensions of these theories.) 



3 Adiabatic Mechanism in an Expanding Phase 

The key element of the adiabatic ekpyrotic mechanism is the rapid evolution of the equation of 
state parameter e; this growth is responsible for the generation of a spectrum of adiabatic modes. 
Interestingly, the scale factor remains nearly constant in the process, and hence plays no essential 
role. As a result, we expect that the adiabatic ekpyrotic mechanism can be generalized to a case 
where the scale factor is slowly increasing. In this Section we show that this is indeed possible. 

Though this may seem like an academic exercise, there are important implications in extending 
the original set-up, especially for embedding this mechanism in broader cosmological scenarios. 
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As reviewed above, in the contracting adiabatic mechanism the transition phase is a dynamical 
attractor [37]. However, the cosmological evolution prior to the onset of the transition phase is 
not necessarily an attractor solution. In fact, we know that at the earliest times when the field lies 
on the potential plateau, the background cosmology must be unstable — any additional kinetic 
energy in the field will blueshift and can lead to kinetic domination. In an initially expanding 
universe, the kinetic energy will instead redshift. In this way, we expect an expanding transition 
phase solution to be stable for all time. 



3.1 The Expanding Transition Solution 

The construction of the solution is quite similar to the contracting case. As before, we may ignore 
the background evolution to first order so the scalar field equation of motion reduces to ([6]), with 
a solution given again by 

log I J-r^r\i\ j . (11) 




where — oo < t < 0. Since the cosmological background is to first order irrelevant for the scalar 
field evolution, it is natural that the approximate solution ([7]) applies irrespective of whether the 
universe is contracting or expanding. In the same way as before, we can integrate the H equation 
to obtain 

H(t) = ^- t +H - (15) 
The Friedmann equation d9|) constrains the magnitude of the integration constant H as before, 



Vo 



Note this time we have chosen the positive-definite quantity H to appear as +H in (15). Hence, 
at sufficiently early times, H(t) ~ Ho is approximately constant and positive, corresponding to an 
expanding de Sitter universe. The expanding transition solution will be similar to the contracting 
case as long as H(t) ~ H , we will have z ~ 1/t 2 , and scale invariant perturbations will be 
generated. However since Hubble is positive during this time, perturbations are generated on a 
slowly expanding background. 



The transition phase is once again followed by an ekpyrotic scaling phase, with H e ^(t) ~ 2/c 2 t. 
During this phase, the Hubble radius shrinks slowly, and the scale invariant modes are gently 
pushed outside the horizon. There is, however, a key difference: in the expanding case, H(t) 
changes sign as the universe evolves from the (expanding) transition phase to the (contracting) 
ekpyrotic scaling phase. 

By definition, the transition phase proceeds as long as the Hubble parameter is nearly constant 
and the Hubble friction term can be neglected in the scalar field evolution. These conditions are 
only satisfied for a finite time. First note that the constant term dominates in the expression for 
H(t) until 

2 

^end— tran ^beg— ek 2 tt ' (1^) 

C SJQ 

which corresponds to the time when H vanishes, and the universe transitions from expansion 
to contraction. As the subscripts indicate, this marks the end of the transition phase and the 
onset of the ekpyrotic scaling phase. Likewise, the transition phase is also finite in the past. 



The solution (14) for cj)(t) neglected gravity, which is a poor approximation for sufficiently large 
positive where the potential is flat and Hubble damping is important. More precisely, it is 
straightforward to show that the approximation Hep <C cVoe~ c ^^ Mpi /Mp\ used in (p| is consistent 



as long as t > tbeg-tran, where 

1 

^beg— tran TT • (1^) 

-"0 

Summarizing, the transition phase solution is valid for 

1 2 

tj ^beg— tran ^ t <i t C nd— tran o TJr ■ (. **) 

Interestingly, the length of the transition phase lasts less than a Hubble time, and thus represents 
a small amount of the total cosmological evolution. 



We can calculate the equation of state parameter e using our expression for H in (15): 

H c 2 /2 c 2 /2 



(20) 



H 2 (1 + H cH/2) 2 (1 - i/t cnd -tran) 2 ' 

At early times, \t\ 3> |4nd-tran|, this expression reduces to e ~ 1/t 2 , characteristic of the transition 
phase. Hence we expect that this will lead to the generation of scale invariant modes. At late 
times, \t\ <C |t e nd-tran|) w e have e w c 2 /2, characteristic of the ekpyrotic scaling phase. 
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Note that the equation of state parameter diverges at t — t e nd-tran — the Hubble parameter 
vanishes at that time while H remains finite. At first sight this may seem worrisome, since the 
momentary divergence of e could potentially disrupt the scale invariance of the perturbations. 
This turns out to be an unfounded fear. We will show in Sec. [4] by numerically integrating the 
perturbation equation that the spectrum for ( is unaffected by the momentary blow-up in the 
equation of state. Intuitively, this is because the time scale over which the equation of state 
diverges is very short, i.e. much less than a Hubble time. 



For completeness, we can integrate (15) to obtain an expression for the scale factor 



a{t) ~ (-tf c2 e Hot ^1 + 4 M-t) + H t + . . . , (21) 

where the Taylor expansion is a good approximation during the transition phase. The scale factor 
is therefore nearly constant during the transition phase, as in the contracting example, and remains 
finite throughout the evolution. 

It is worth emphasizing that while the background is expanding during the transition phase, this 
is decidedly not an inflationary scenario. During the expanding transition phase, the equation of 
state evolves rapidly in time (e(i) ~ l/^ 2 ), while the scale factor is nearly constant (a(t) ~ 1), in 
stark contrast to slow-roll inflation. 



3.2 Curvature Perturbations in the Expanding Phase 



We are now in a position to check that a slowly-expanding transition phase leads to a scale invariant 
spectrum of density perturbations. We begin with an analytic calculation, before verifying the 
results numerically in Sec.|4j For this purpose, we work in comoving gauge, where the spatial slices 
are constant density (5<p = 0) hypersurfaces, and the spatial metric is given by hy = a 2 (t)e 2 %. 
In this gauge, the variable ( represents the curvature perturbation on spatial slices. Its action at 
quadratic order is given by 



So = - d xdr z 
2 J 
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vc 



(22) 



where z = ay2e as before, and primes denote derivatives with respect to conformal time r. The 
resulting equation of motion for the Fourier modes C,k is 



C + 2-G + rc fc = o. (23) 

It is convenient to instead work in terms of the canonically normalized variable v = z ■ (: 

< +(k 2 --)v k = 0. (24) 



z 



If z" jz ~ 2/r 2 , as in inflation, then this equation is well-known to yield a scale-invariant spectrum 
for ( at long wavelengths. Furthermore, in this case the growing mode of ( goes to a constant 
as k — > 0. In this limit, ( may be interpreted as a homogeneous perturbation to the scale factor 
and may be locally absorbed by a spatial diffeomorphism [35]. The background is therefore a 
dynamical attractor. 



In our case, it follows from (20) and (21) that 

z(t) = (-t) 2/c2 e^*— . (25) 

-L "I ''end— tran 

Thus z inherits the singular behavior of e at t — t e nd-tranj though as mentioned earlier this will 
not pose a problem. One way to see this is to note that t = t e nd-tran is a regular singular point 



of (|23j). 

During the transition phase, the scale factor is nearly constant, hence conformal time and cosmo- 



logical time are approximately the same: t w r. Deep in the transition phase, \t\ ^> |£ e nd-tran|, (25) 
therefore implies 

— ~ 2 _ <!_ 

Z (t ^end— tran) T 

Assuming the usual adiabatic vacuum, the solution is 



which gives a scale invariant power spectrum for the curvature perturbation, ( = v/z, in the long 
wavelength limit (k\r\ 1): 

p < = J_ ^| Cfc |2 = ^0 ( } 

C 2tt 2 ISfe| 48vr 2 V ; 
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Moreover, in this limit £ approaches a constant, indicating that the expanding transition phase is 
also an attractor. 



Since modes freeze out when k\r\ ~ k\t\ ~ 1, the range of scale invariant modes is set by the 



duration of the transition phase. From (17) and (18), we deduce that 



^be 



k ■ 



«end- 



c 
2~ 



(29) 



Equations (28) and (29) agree exactly with the contracting case EZ] • Of course, this is not 
surprising, since away from the singular point t = t e nd-tran, z does not know the difference between 
expanding and contracting solutions. We now check such claims through numerical integration. 



4 Numerical Verification of Scale Invariance 

In this Section we check the analytical results derived above by numerically integrating the evo- 
lution of (. In the process we will reassure ourselves that the momentary singularity in e does not 
spoil the scale invariance of the modes created during the transition phase. 



Our starting point is the evolution equation (23) for the mode functions This equation is 



better behaved than that for Vk, since, as mentioned earlier, the singularity at t = t e nd-tran is a 
regular singular point in this case. As initial conditions we impose the adiabatic vacuum choice 

4(Ti) - • (30) 

with Ti chosen so that kWA 3> 1 for all modes of interest. 



To proceed, we need to substitute an expression for z(r) into (23). We will do this in two ways. 



First, following a quasi-analytic approach, in Sec. 4.1 we use the analytic form given by (25). This 



expression is of course not exact, since various approximations went into deriving it. Second, in 



Sec. 4.2 we will redo the calculation more precisely using an exact form for z(r), which will be 
obtained by numerically integrating the background equations of motion. 
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4.1 Integration with Analytic z{r) 



Let us start with the quasi- analytic approach, using the approximate form for z{r) given by (25). 



As it stands, (25) gives z as a function of cosmological time t, whereas we need it in terms of 
conformal time r. Note that the factors of c/(l — £/t nd-tran) and (— t) 2 ^ 2 are important in the 
transition phase and the ensuing contracting ekpyrotic phase, respectively. During both of these 
phases, cosmological time and conformal time are nearly the same, so a good approximation is 
to replace t with r in these factors. The factor of e Hot , however, is important at early times and 
makes the scale factor differ from unity — conformal time and proper time are therefore much 
different in this case. Approximating a(t) ~ e Hot at early times, we can integrate dr = dt/a(t) to 
obtain 

Thus we obtain the following analytic expression for z(r) 

c(-t\ 2 I (? 
(1 - Hot) (1 - r/W-traa) 

Using this analytic expression for z(r) we are now in a position to integrate the equation of motion 
for (. In fact, we will do the integration for both the expanding and the contracting case, where 



the contracting case corresponds to letting H — > —H in (32). In both cases, we take c = 200 



H = 5 x 10~ 4 , and integrate (23) over the range of modes .02H < k < 2 x 10 A H and over the 
time interval —.hH^ 1 < t < —5 x 10~ 10 H ~ 1 . These parameters are chosen so that the shortest 
wavelength modes will have just left the Hubble horizon by the end of the evolution, while the 
background solution will be firmly within the scaling ekpyrotic phase. 

The resulting power spectra are shown in Fig. [2] We find that there is a range of scale-invariant 
modes spanning roughly 4 decades in fc-space, which is in order-of-magnitude agreement with our 



analytical result (29). Note that the range of scale invariant modes is slightly shorter than that 
found in the numerical analysis of [35] , which is due to our including the extra factor 1 /(l — Hqt) in 
z{t). We see that the power spectrum for ( is indistinguishable in each case, which confirms that 
the adiabatic ekpyrotic generation mechanism works equally well on an expanding background. In 
particular, this confirms that the divergence in e at t — t e nd-tran, corresponding to the transition 
from expansion to contraction, has no effect on the perturbation spectrum. 
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Figure 2: Numerical computation of the power spectrum k 3 ^ 2 (k vs. k in both the expanding and 
contracting transition phase scenarios using an analytic expression for z(t). This confirms that 
the spectrum is insensitive to whether the transition phase is contracting or expanding. 

4.2 Numerical Solution of the Full Background Equations 



The above analysis relied on certain approximations in deriving (32). We neglected Hubble friction 



in the evolution for </>, but this must eventually break down at early times when the background 



is approximately expanding de Sitter space. To check that such corrections to (32) do not spoil 



scale invariance, in this Section we numerically solve for the background evolution to obtain an 
exact result for z{t). 

This can be done using the Hamilton- Jacobi formalism. Using the chain rule, the H equation 
implies 

H* = -J^, (33) 



hence the Friedmann equation can be rewritten as 

1 



SM^H 2 = + V{4>) = 2M^H% + V{4>) 



(34) 



In this way, is thought of as the clock tracking the background evolution. We first numerically 



solve this equation to find H(<p), and in turn obtain <p(t) by integrating (33). At the end of 
the day, this gives us the Hubble parameter as a function of time H(4>(t)), from which we can 
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Figure 3: Comparison of the power spectra for the exact numerical calculation and the integration 
using our analytic expression for z(t). The curves are in excellent agreement. 



extract a(t). In practice, we perform the integration over a sufficiently broad range of field values 



and correspondingly large time interval. Specifically, we have solved (34) over the field range 
0.07 Mpi <<p< -M Ph setting #(0.07 Af P i) = H . Meanwhile, in integrating Q to obtain <p(t) 
we fix the integration constant corresponding to a shift in time by demanding that matches the 



analytic solution (14) 



Substituting everything into z = ay2e, we can numerically integrate (23) to obtain the power 
spectrum. The result is shown in Fig. [3^. For comparison, Fig. |3)d shows the quasi-analytic result. 
Computational constraints forced us to use slightly different integration parameters from those 
used in Sec. 4.1 Figure^ was obtained using c = 200 and H = 10~ 3 , with integration ranging 



over 



-.5#n 



< r < 



TO 8 H 1 and .01H < k < 10 4 H . We see that the plots show very good 



agreement. 



5 Starting From Rest 

The analysis of Sec. [3] assumed a scalar field starting from rest (0 = 0) in the asymptotic past. 
In this Section we explore a broader range of initial conditions, corresponding to non-zero initial 
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kinetic energy for the scalar field, and the impact on the range of scale invariant modes. Specifi- 
cally, our goal is to derive an expression for the length of time spent in the transition phase (and 
thus the range of scale invariant /c-modes generated) as a function of the initial energy. 



We start with the equation of motion for the scalar field, assuming as before that the Hubble 
damping term is negligible 

4> + « . (35) 



This equation admits a first integral of motion 

1 ■ 



2 + V(<t>) = E 



(36) 



where E is of course the total energy of the field. In Sec. [3j as well as in earlier work EZ] , E 
was taken to be equal to Vq at the onset of the transition phase, while remaining agnostic about 
the prior evolution. Here we want to consider the more general case E < Vq. In other words, E 
corresponds to the value of the potential from which the scalar field starts at rest. 



Equation (36) implies 



= -v/2(£-K(0)) , 



(37) 



where we have chosen the negative branch of the square root, corresponding to the field rolling 
downhill. Substituting our potential Q, this can be integrated explicitly 



- t 



C(f>/2M P 



^ arcsin ( J^^e^" 



(38) 



y/2 [V + (E- V )e^ M ^] c^2(V - E) " \ V \ o 

where the last step assumes e c4> / 2 < y/V /(V -E). This can be inverted to obtain a solution for 
(pit), valid as long as <p < (2/c) log y/V /(V - E): 



2Mi 



pi 



log 



Vn 



V n -E 



V -E . . 
' 1 ' ~2M| 1 ~ C 



(39) 



This expression generalizes (14) to span a range of initial field energy, or equivalently initial 



position on the potential where the field is released from rest. As a check, this reduces to (14) 
in the limit E — > Vq. Moreover, at late times, when \t\ becomes sufficiently small, (39) also 



approaches (14), confirming that this is an attractor. 
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It is straightforward to calculate H 



Vn-E 




H =-2Mf* --^T COt 'UHri^l'll " ( 4 °) 

Similarly, substituting the above solution for cf>{t) in the Friedmann equation, 3H 2 Mp { = + 
V(<j)) = E, with V(4>) = Vo(l — e~ c ^), we obtain H 2 = E/3Mp V The equation of state parameter 
is therefore given by 

H 3(V - E) ( j v -E , 

e =-jp = — E— cot v^ir c|11 ■ ( } 



This generalizes (20 ) to a broader range of initial conditions, and matches (20 ) in the limit E — > Vq, 



as well as at late times \t\ —> 0. 

Clearly the range of scale invariant modes will depend on the initial conditions, and we can 
already expect a more restricted range as we move away from the case E = Vq studied earlier. 
The relevant quantity to assess the shape of the power spectrum is the time-dependent mass term, 



z"/z, appearing in the mode function equation (24). Since a(t) is nearly constant during the phase 
of interest, we have 

z" 1 d 2 Je c 2 (V -E) 2 . , .„ ... . . . 

- V^ csc ' ( \I^II^c\t\ I • ( 42 ) 



z dt 2 Mj 



pi 




Modes will be scale invariant provided they freeze out when z"/z ~ 2/t 2 . Recalling the Taylor 
expansion esc 2 x ~ x~ 2 + 1/3, this will be the case whenever c(-t)y/(y - E)/2 < 1. This 
condition is approximately satisfied for t < tbeg-tran, where 



the ^ n ^ ~ly Vo^e ■ (43) 

In other words, tbeg-tran marks the onset of scale invariant mode production. Although this 
expression naively diverges for E = Vq, we of course can only make ibeg-tran as large as tbeg-tran = 
— 1/Hq — at earlier times, the fast-roll approximation assumed here breaks down. More carefully, 
we have 



tbeg-tran - max ( ~\^y^, ~ | ■ (44) 
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The scale invariant phase has finite duration. Much like the transition phase studied earlier, it 



comes to an end when the approximation H ~ v/ 'E /3M|j breaks down. Integrating (40) in the 
limit c(-t)y/(V - E)/2 < 1 gives 



H(t) 



E 2 
3M^ + dH + "' 



(45) 



where the ellipses include terms that become increasingly small as t — > 0. It follows that the scale 
invariant or transition phase concludes at 



^end- 



E 



(46) 



Combining (43) and (46), the range of scale invariant modes is thus given by 

V 



t 



beg— tran *- 




(47) 

'■end— tran 

Here we have assumed the earliest possible onset of the transition phase. Comparing to our 

c 2 /2, we see that solutions with E < Vq lead to a narrower 



previous answer (29), k max /k 



range of scale invariant modes, as expected. Note that (47) agrees with (29) in the limit E — > Vq, 
as it should. 



6 Phase Space Analysis 



In Sec. 3.2, as well as in earlier work [351 [37], the adiabatic ekpyrotic evolution has been argued to 
be a dynamical attractor because the growing mode for ( goes to a constant in the long wavelength 
limit |45j. This argument, while true, only applies at the perturbative level, i.e. for sufficiently 
small deviations from the background solution. In this Section, we study the issue of stability 
more broadly, by performing a phase space analysis for a wide range of initial conditions. This 
will allow us to determine the breadth of the basin of attractor both in the expanding and in the 
contracting branches of the transition phase. 

For this purpose, we consider curves in the (0, <fi) phase plane parameterized by iV = log a. The 
Friedmann constraint and the scalar field equation of motion imply the autonomous system: 

d0 ^ 
diV ~ H' 
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where the Hubble parameter is understood as a function of <f> and 4>: 

The choice of sign for the square root corresponds to the choice of expanding or contracting branch. 
Cosmological solutions are given by the integral curves of the vector field (d0/diV, d0/diV). 



6.1 Definition of Attractor 

Let us be precise about what we mean by attractor behavior. As curves approach each other in 
the ((f), cj)) plane, their respective values of <p(N) and (f>(N) get closer and closer together, up to a 
constant relative shift in N. We have the gauge freedom to fix the initial values of N such that the 
solutions coincide. In this sense, an attractor solution is identified by neighboring curves flowing 
towards it. We will see that this is the case for the transition phase and the subsequent ekpyrotic 
scaling phase. 

It is worth pointing out that this focusing of trajectories occurs in the (0, 0) plane, but not in 
the physical phase space. Minisuperspace models are Hamiltonian systems, and consequently 
Liouville's theorem forbids reduction of phase space volume. The resolution of this apparent 
discrepancy is of course that = a 3 0, rather than <ft, is the momentum conjugate to <ft i n the 
canonical formalism. 

This can be seen more generally by considering a 2n dimensional phase space. Minisuperspace is 
a constrained Hamiltonian system, with the Friedmann equation defining a 2n — 1 dimensional 
constraint hypersurface on which the Hamiltonian vanishes. This hypersurface is foliated by the 
gauge orbits corresponding to the Hamiltonian flow. The 2n — 2 dimensional hypersurface which 
is transverse to this Hamiltonian flow is the space of classical trajectories; this is the symplectic 
reduction or Marsden-Weinstein quotient M = "H _1 (0)/1R, where 'H^ 1 (0) is the locus where the 
Hamiltonian vanishes [57J EE] • There is a well-defined pullback of the symplectic form in the total 
space to this reduced phase space. Darboux's theorem tells us that there are local coordinates 
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(pi, q l ) such that the symplectic 2-form is 

n 

oj = dpi A dg* . (50) 

i=i 

For a model with action Q in the ADM formalism we can choose coordinates to write this 
symplectic form as 

n-l 

u = m a dt + d ^ A • ( 51 ) 

i=l 

In this language, the Friedmann equation is the constraint that the Hamiltonian vanishes 1-1 = 0. 
The total phase space carries a natural invariant measure constructed from u. We may pull back 
this symplectic form to the reduced phase space M and construct a natural measure fl This is 
the Gibbons Hawking Stewart measure [67]. 



Note that on the reduced phase space, this measure is conserved. Moreover, in the particular 
case of a spatially-flat (k = 0) universe we can write the equations of motion in such a way that 



they are independent of a, as we did in the single-field autonomous system (48). Focusing on this 
single-field case for concreteness, we may then choose to parameterize the reduced phase space by 
and (f> [HH]- Choosing to treat these as Euclidean coordinates introduces the measure 

fi = d0 A d0 , (52) 

but this of course need not be conserved by the evolution. As mentioned earlier, this non- 
conservation can be traced to the fact that is not the momentum canonically conjugate to 



. Using (52) is nevertheless a sensible thing to do — as mentioned above, we are interested in 



trajectories with the same values of and 



6.2 Numerical Results 

Figure [4] shows the vector field (d0/diV, d0/diV) given by (48) for the expanding and contracting 
cases, respectively, along with some numerically integrated curves. The black curve in each case 
corresponds to the analytic solution, with the transition phase occurring between the dotted lines. 
The chosen parameters are c = 100 and Vq = 10~ 4 . Figure [5] zooms in on the transition phase. 
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(a) Expanding (b) Contracting 

Figure 4: Phase portrait for the expanding (a) and contracting (b) cases, for c = 100 and Vq = 
10~ 4 . The analytic solution is denoted by the black curve, with the transition phase taking 
place between the dotted lines. Colored dashed lines are particular numerical solutions to the 



system (48). Note that the arrows point in the direction of increasing time. This Figure confirms 
that the analytic solution is an attractor for a variety of initial conditions in both cases. However, 
at large positive where V(4>) ~ Vo, the expanding solution is also an attractor, while the 
contracting solution is a repellor, due to the asymmetry between expanding vs. contracting de 
Sitter space. 
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(a) Expanding (b) Contracting 

Figure 5: Detailed view of Fig. |4j zooming in on the transition phase. Both contracting and 
expanding transition solutions are attractors, but the expanding case has a slightly larger basin 
of attraction. 
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As is clear from Fig. [5j the transition solution is an attractor, both in the expanding and contract- 
ing cases. This confirms our earlier perturbative claims based on the long-wavelength behavior 
of (. The expanding solution has a slightly larger basin of attraction, as expected. In particular, 
note that in the contracting case curves with (ft sufficiently greater than zero are forced away from 
the transition solution while in the expanding case, such curves are driven toward the slow-roll 
solution and then follow the transition solution. 

More globally, Fig. [4] shows that in the contracting case curves are repelled from the 'slow-roll' 
solution at large <fi. This is due to the susceptibility of contracting de Sitter to kinetic domination. 
In the expanding case however, we see the opposite behavior — trajectories with significant initial 
kinetic energy are driven toward the slow-roll de Sitter solution before undergoing the transition 
phase evolution. To summarize, both the expanding and contracting transition solutions are 
attractors for some range of initial conditions, but the expanding solution has a larger basin of 
attraction. 



7 Fusion of Slow-Roll and Transition Solutions 

One of our motivations for considering an initially expanding cosmology was the instability of the 
contracting solution at sufficiently early times. The adiabatic ekpyrotic mechanism requires that 
the transition phase starts out with a sufficiently small equation of state parameter, e< 1, which 
implies that prior evolution will generically be unstable to kinetic domination in the contracting 
case. 

As shown in Sec. |6j this is true, in particular, for the simplest potential ([I]) — the pre-transition 
phase evolution corresponds to a contracting de Sitter universe, which is clearly unstable. The 
expanding branch, on the other hand, is better behaved, as it extrapolates backwards in time to 
an expanding de Sitter space. Thus, in choosing an initially expanding universe, we have greatly 
expanded the basis of attraction. (It is worth emphasizing that even in the initially expanding 
case, the universe is driven towards a contracting ekpyrotic scaling phase.) 

As mentioned in the Introduction, our fiducial potential ([T| serves as the simplest illustration of 
our mechanism. More generally, this form need only hold approximately during the transition 
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phase, corresponding to a small range in field space, A0 <C Mpi, and there is ample freedom in 
choosing the global form of the potential. That said, if we take the potential Q at face value, then 
at sufficiently early times the evolution in the expanding case should correspond to an inflating 
space-time. In this Section, we take this possibility seriously and study the transition between the 
initial de Sitter phase to the transition phase. We will see that <p(t) evolves smoothly between the 
two regimes. 

To see this, let us split the evolution into two regimes: 



I The 'transition' regime, < 0t> where Hubble friction is negligible and the scalar field 
equation of motion reduces to 

4> + V^ = 0. (53) 
II The 'slow-roll' regime, <fi > 0t> where the equation of motion is approximately given by 

3H<j) + = , (54) 

with H ~H = y/Vo/3^. 



The field value delineating these two regions, 0t, can be estimated as the value of at the onset 



of the transition phase. Combining (14), (16) and (18), we obtain 




4>T = 0(ibeg-tran) = log W -C ~ log C , (55) 



C 

where in the last step we have used c^> 1. 

We are now in a position to argue that the approximate solutions in regions I and II match each 
other smoothly at <p = 0t- Since the potential is certainly continuous, this amounts to showing 
that the kinetic energy also matches continuously. The kinetic energy of the transition solution is 



easy to write down, as we have done it several times already, and is just the first integral of (53): 

^Ln « V e-^ . (56) 



Likewise, the kinetic energy in the slow-roll regime follows from (54): 

2 Vslow-roll~ 
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Lw-ron - - V c>e-^ . (57) 



It immediately follows that (56) and (57) are equal at = fa, as we wanted to show. This 
result tells us that trajectories starting in the slow roll region, <fi > fa, are quickly driven to the 
attractor slow-roll solution and reach fa with precisely the correct kinetic energy as assumed by 
the transition solution. (In the notation of Sec. |5j solutions that emerge from the slow regime 
generically reach the transition phase with E ~ Vo.) 

Note that the spectrum of fluctuations also matches smoothly between the two regimes. The 
power spectrum of modes generated during the inflationary epoch is given by 

P m - _i ~ <^L (58) 



which agrees up to an order unity factor with ( 28 ) . 



8 Conclusions 

The adiabatic ekpyrotic mechanism is the unique, non-inflationary single-field mechanism with 
unit sound speed that generates a scale invariant spectrum for ( on an attractor background. 
While originally proposed assuming a contracting universe, in this paper we have shown that the 
mechanism works equally well on an initially expanding background. The evolution consists of an 
"expanding transition phase" , followed by a contracting ekpyrotic scaling phase. 

We have shown, both through analytical arguments and exact numerical integration, that the 
power spectrum for ( is indistinguishable from the contracting version of the mechanism. This 
confirms that the perturbation spectrum is to a good approximation insensitive to whether the 
background is expanding or contracting during mode production — scale invariance relies on a 
rapidly- varying equation of state parameter while the scale factor is nearly constant. 

By performing a phase space analysis for both the expanding and contracting branches, we have 
verified that the transition phase and subsequent ekpyrotic scaling phase are attractors. The 
basin of attraction is broader in the expanding case, since intuitively any additional kinetic energy 
present at early times gets redshifted, instead of blueshifted, in this case. 



For the simplest potentials considered here, the evolution is an expanding de Sitter space-time 
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asymptotically in the past. Taking this precursor inflating phase seriously, we have shown that 
the transition from inflation to the adiabatic ekpyrotic phase is smooth, with the scale invariant 
spectra generated in each phase matching at the transition, up to order unity coefficient. 

As mentioned in the Introduction, the degeneracy with inflation is broken at the three-point level. 
Unlike the nearly Gaussian spectrum of inflation, the rapidly- varying equation of state e ~ 1/r 2 
characteristic of the adiabatic ekpyrotic phase leads to large non-Gaussianities on small scales. For 
the simplest lifted exponential potential 0, this results in a breakdown of perturbation theory, 
both at the classical and quantum mechanical levels. As shown in detail in [37] in the contracting 
case, however, this perturbative breakdown can be avoided for more general potentials of the 
form ([3]), but this comes at a cost — the range of scale invariant perturbations is now limited, 
spanning a factor of 10 5 in k space. These considerations should carry over to the expanding case. 
We leave a careful exploration of this issue to future work. 

The expanding version of adiabatic ekpyrosis presented here greatly expands the realm of larger 
cosmological scenarios in which this mechanism can be embedded. Because the basin of attraction 
is broader compared to its contracting counterpart, the expanding mechanism is less sensitive to 
the details of the prior evolution. More importantly for model-building, the universe is expanding 
during mode production and is contracting later on, which naturally suggests embedding our 
mechanism in a cyclic scenario [H] . It is conceivable that tying our mechanism to the present phase 
of cosmic acceleration may also explain why the 10 5 scale invariant modes fall in the observable 
window. We are currently investigating this possibility. 
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